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CSE&A  fCTHODOLOGY  DEVELOPMENTS 


I 

INTRODUCTION 

The  Military  Effectiveness  Office  {DTNSRDC  1806)  has  been  support¬ 
ing  the  NAVSEA  (PMS  393)  SSN  Combat  Systems  Engineering  and  Analysis  (CSE&A) 
Program  by  producing  the  Attack  Submarine  Development  Plan  (ASDP).  As  part 
of  this  effort,  a  methodology  for  evaluating  the  military  payoffs  of  adding 
various  advanced  equipments  to  the  baseline  SSN  716  has  been  developed.  This 
effort  Is  directed  at  Improving  and  extending  the  SSN  CSE&A  evaluation  models 
dealing  with  three  areas  of  SSN  performance  Including: 

1)  SSN  effectiveness  In  a  moving  barrier  mission 

2)  SSN  effectiveness  In  an  area  clearance  mission 

3)  SSN  effectiveness  In  multiple  engagement  scenarios. 

This  report  documents  the  results  of  the  effort  up  to  the  time  when 
the  task  was  terminated  due  to  funding  cutbacks.  The  barrier  model  was  com¬ 
plete  in  Its  development  except  for  the  testing  of  the  model  which  was 
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ongoing.  The  area  clearance  algorithm  was  developed  but  the  model  necessary 
for  generation  of  results  was  not  completed.  No  work  was  performed  on  multi¬ 
ple  engagement  effectiveness^ 

BARRIER  MODEL 

The  barrier  model  was  used  to  compute  the  probability  that  a  sub¬ 
marine  transiting  a  channel  of  width  Lg  at  speed  Vj  will  be  detected  at 
some  time  during  Its  passage  by  a  barrier  submarine  which  patrols  back  and 
forth  across  all  or  part  of  the  channel  (Lp)  at  speed  Vg.  The  transltor 
Is  assumed  to  be  travelling  perpendicular  to  the  barrier  since  this  Is  his 
best  tactic  to  minimize  detection  opportunity. 

For  initial  positions  of  transltor  and  barrier  submarine  shown  In 
Figure  1A,  the  relative  track  of  the  transltor  In  barrier  space  (l.e.,  the 
barrier  submarine  Is  the  fixed  center  of  this  moving  coordinate  space)  is 
shown  In  Figure  IB.  The  Initial  direction  of  the  barrier  submarine  Is  Indi¬ 
cated  by  the  arrow  above  the  B  which  denotes  the  barrier  submarine. 

This  relative  path  is  only  one  of  an  Infinite  number  of  possible 
paths  related  to  the  Initial  positions  of  the  transltor  and  barrier  subma¬ 
rine.  The  transltor  submarines  Initial  position  Is  uniformly  distributed  at 
some  maximum  range  across  the  total  barrier  length.  The  barrier  submarines 
Initial  position  Is  uniformly  distributed  over  the  barrier  patrol  length  with 
Initial  direction  of  travel  either  to  the  right  or  left.  This  Is  true  except 
at  the  end  points  of  the  patrol  at  which  only  one  direction  Is  possible. 

The  barrier  effectiveness  then  is  the  mean  probability  of  detection 
over  these  Infinite  combinations  of  paths. 

The  probability  of  detection  along  a  particular  path  Is  computed 
using  a  "relaxation  time"  model  In  order  to  deal  with  the  correlation  between 
samples  taken  continuously  along  the  path.  As  stated  in  Reference  1,  let 
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Figure  1A  INITIAL  TRANSITOR  AND  BARRIER  SUBMARINE  POSITIONS 
AND  BARRIER  DIMENSIONS 


Figure  IB  RELATIVE  TRACK  OF  TRANSITOR  IN  BARRIER  SUBMARINE  SPACE 
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x(t)  be  a  0-mean  stochastic  process  and  define  E(t)  *e(t)  -  x{t)  to  be  the 
time  dependent  signal  excess  where  e{t)  Is  the  "signal  excess"  computed  from 
the  relative  position  of  target  and  receiver  and  by  the  FOM  of  the  sensor. 
The  probability  of  detection  Is  then  a  first  passage  probability  for  E(t) 
through  the  level  0.  The  simplest  case  Is  where  x(t)  is  constant  over  Inter¬ 
vals,  being  chosen  Independently  from  a  normal  distribution  with  mean  0  and 
variance  o2  in  each  Interval.  If  the  Intervals  all  have  the  same  length,  the 
model  is  called  the  "relaxation  time"  model.  In  this  case  the  probability 
that  no  detection  occurs  in  the  1th  interval  Is 

q|  =  1  -  ♦ (e(t| )/o ) 

where  by  definition  the  maximum  of  e(t)  in  the  1th  interval  occurs  at  t|  and 
♦  Is  the  cumulative  normal  function.  The  probability  of  detection  along  any 
path  Is  then 


p  =  1  -  niq1 

where  the  product  extends  over  whatever  Intervals  are  involved  In  the  time 
period  necessary  to  complete  one  relative  track. 

Having  to  break  up  the  track  Into  equal  segments  for  the  "relaxa¬ 
tion  time”  model,  also  allows  a  reduction  In  the  number  of  paths  needed  to  be 
considered  and  an  Increase  In  computational  efficiency  when  the  length  of  the 
segments  Is  chosen  carefully. 

The  barrier  length  Is  divided  Into  N  equal  segments  that  roughly 
correspond  In  length  to  the  distance  travelled  by  the  barrier  submarine  In 
one  relaxation  time  period.  The  relaxation  time  period  Is  related  to  the 
macro-environmental  fluctuations  of  the  medium.  The  barrier  submarine  re¬ 
quires  a  time  tg  to  travel  a  distance  equal  to  the  length  of  one  of  these 
segments.  The  transltor  meanwhile  travels  a  distance  V-j-tg  In  the  time 
tg. 
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The  transltor  starts  his  run  at  a  distance  (Rmax)  from  the  bar¬ 
rier  outside  the  maximum  possible  range  of  detection  (that  range  for  which 
the  detection  probability  Is  at  or  near  zero).  The  transltor  travels  M  range 
segments  of  length  (V^tg)  to  reach  the  barrier  line  from  this  initial 
range.  After  the  transltor  has  moved  from  Rmax  to  the  barrier  line,  the 
barrier  submarine  has  patrol ed  M  barrier  segments  back  and  forth  from  Its 
Initial  position.  At  the  barrier  line  the  barrier  submarine  and  the  transl¬ 
tor  will  be  separated  by  some  number  of  barrier  segments,  as  shown  In  Figure 
IB. 


In  barrier  space  we  have  developed  a  grid  of  points,  2N+1  by  2M 
(Rmax  to  -  Rmax)  where  the  relative  transltor  track  can  be  reconstructed 
by  appropriately  moving  diagonally  from  one  grid  point  to  the  next  and  where 
changes  In  the  direction  of  diagonal  movement  correspond  to  the  reversal  of 
barrier  submarines  direction  at  the  end  of  the  patrol  length.  (See  Fig¬ 
ure  IB.) 


If  Np  segments  out  of  the  total  N  barrier  segments  are  partolled 
(N  -  Np  must  be  a  multiple  of  2  to  ensure  that  the  patrol  length  Is  cen¬ 
tered  in  the  barrier)  then  there  are  Np  +  1  endpoints  of  the  patrol  seg¬ 
ments.  The  barrier  submarines  Initial  position  Is  allowed  to  be  at  any  of 
these  points.  Initially  travelling  In  either  direction,  except  at  the  end 
points  of  the  patrol  length  where  only  one  direction  Is  possible.  There  are 
therefore  2Np  Initial  starting  points  (taking  into  account  position  and 
direction)  that  are  possible  for  the  barrier  submarine.  The  transltor  sub¬ 
marine  may  start  at  any  one  of  the  N  +  1  points  of  the  barrier  length  from  a 
distance  of  M  range  segments  from  the  barrier  line.  The  total  number  of 
paths  to  be  considered  Is  2Np(N-<-l),  accounting  for  all  possible  combina¬ 
tions  of  Initial  positions  for  transltor  and  barrier  submarines. 

The  barrier  model  permits  utilization  of  aspect  dependent  source 
level  data.  In  barrier  space,  the  transltor  progresses  along  the  relative 
track  with  heading  perpendicular  to  the  barrier  line  (Figure  IB).  Thus  for 
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any  point  In  the  grid,  the  aspect  angle  Is  the  aspect  presented  to  the  center 
of  the  grid  (barrier  submarine).  The  probability  of  detecting  the  transitor 
when  moving  diagonally  from  one  grid  point  to  the  next  is  the  highest  proba¬ 
bility  of  detection  (i.e.,  highest  signal  excess)  in  this  segment.  In  com¬ 
puting  detection  probabilities  with  aspect  dependent  source  level  data,  all 
terms  of  the  signal  excess  equation  remain  constant  except  for  the  source 
level  and  propagation  loss  terms.  The  highest  signal  excess  in  a  segment  of 
relative  track  is  found  by  determining  the  maximum  of  source  level  minus 
propagation  loss  for  all  points  in  that  segment.  If  source  level  is  assumed 
independent  of  aspect,  the  maximum  signal  excess  for  a  segment  of  relative 
track  occurs  when  the  propagation  loss  is  at  its  minimum  in  this  segment 
since  propagation  loss  is  the  only  variable  in  the  signal  excess  equation. 

The  probability  of  not  detecting  the  transitor  is  computed  for 
every  grid  point  for  the  diagonal  segments  entering  that  grid  point  from  the 
right  and/or  left,  depending  upon  the  physical  possibility.  Once  the  proba¬ 
bility  of  not  detecting  is  calculated  for  all  grid  points,  an  indexing  scheme 
is  used  to  calculate  all  possible  tracks.  This  scheme  is  presented  in 
Appendix  A. 


When  reconstructing  the  2Np(N+l)  possible  tracks,  most  of  these 
diagonal  segments  are  included  in  more  than  one  track.  The  computational 
efficiency  of  the  model  is  therefore  improved  since  the  probability  of  not 
detecting  in  any  segment  is  only  calculated  once  even  though  it  may  be  used 
in  numerous  track  reconstructions. 

The  overall  probability  of  detecting  the  transitor  along  any  rela¬ 
tive  track  is  therefore  l-fyp  (not  detecting  in  segment  1)  and  the  barrier 
effectiveness  is  the  average  probability  of  detecting  the  transitor  for  all 
2Np(N+l)  tracks. 
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BARRIER  MODEL  INPUTS 

Inputs  to  the  model  are  as  follows: 

Barrier  Length  (nm) 

Number  of  Barrier  Segments 
Number  of  Patrol  Segments 
Velocity  of  transltor  (knots) 

Velocity  of  barrier  submarine  (knots) 

Background  Noise  received  by  barrier  submarine  sensor  (dB) 
Recognition  Differential  of  barrier  submarine  processor  1  ) 

Source  Level  File  (If  aspect  dependent  data  is  used) 

Propagation  Loss  File 
Choice  of: 

a)  Transltor  submarine  travels  from  R^ax  to  the  barrier 
line  (o)  or 

b)  Transltor  submarine  travels  from  R^x  to  -Rmax 
NOTE:  R^ax  Is  the  maximum  possible  range  at  which 

barrier  submarine  can  detect  transltor 

(probability  of  detection  =*.001). 

BARRIER  MODEL  BENEFITS 

The  barrier  model  as  developed  has  numerous  advantages  over  the  use 
of  the  standard  barrier  probability  of  detection  equations  In  which  the  bar¬ 
rier  submarine  patrols  the  barrier  length  to  within  a  fixed  distance  (W/2)of 
the  barrier  end  points  (W/2  Is  one  half  the  sweepwldth  of  the  sensor).  These 
advantages  Include: 

1)  Use  of  actual  probability  of  detection  versus  range 
and  aspect  rather  than  the  "cookie  cutter"  W/2.  The 
use  of  a  "cookie  cutter”  detection  pattern  consis¬ 
tently  overestimates  the  barrier  probability  of  de¬ 
tection.  In  limited  exercise  of  this  model  It  was 
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seen  that  the  maximum  barrier  effectiveness  occurred 
when  the  patrol  width  was  approximately  one-third  of 
the  barrier. 

2)  Inclusion  of  decorrelation  of  environmental  fluctua¬ 
tions  through  using  the  "relaxation  time"  model  of 
detection. 

3)  Incorporation  of  aspect  dependent  source  level  data, 
if  desired. 

4)  Ability  to  use  any  propagation  loss  function. 

In  addition,  besides  the  computational  efficiency  discussed  earler, 
the  barrier  model  may  be  rerun  using  a  smaller  number  of  barrier  and/or  pa¬ 
trol  segments  without  recomputing  grid  point  probabilities.  This  is  accom¬ 
plished  by  changing  the  range  of  indices  in  the  model  over  which  the  product 
of  terms  Hiqj  is  computed. 

AREA  SEARCH  ALGORITHM 

The  Measure  of  Effectiveness  (MOE)  currently  used  in  the  CSE&A  pro¬ 
gram  for  area  search  is  the  probability  of  detecting  a  target  randomly  situ¬ 
ated  in  area  (A)  given  that  the  searcher  has  a  particular  sweep  width  (W) 
against  the  target  and  travels  at  velocity  (Y)  searching  during  a  time  period 
(T).  The  equation  for  the  probability  of  detection  is: 

PD  =  1  -  exp(-YWT/A) 

The  MOE  to  be  developed  In  this  section  is  the  probability  that 
there  are  no  undetected  targets  In  area  A  after  a  search  of  duration  T.  This 
Is  felt  to  be  a  more  realistic  MOE  for  area  search  than  that  used  in  the 
CSE&A  program. 
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The  targets  for  this  algorithm  are  assumed  to  be  located  randomly 
in  a  two-dimensional  field  with  the  average  number  of  targets  per  unit  area 
equal  to  y.  That  is,  every  small  area,  a,  contains  a  target  with  probability 
ya.  The  probability  of  more  than  one  target  in  a  Is  negligible  if  a  is  suf¬ 
ficiently  small.  Also,  the  number  of  targets  in  non-overlapping  regions  are 
Independent  random  variables.  These  assumptions  are  consistent  with  a  Pois¬ 
son  distribution. 

Let  A  be  the  area  of  a  region  within  the  field  that  is  not  neces- 
sarly  small  and  let  be  the  number  of  targets  in  it.  As  developed  in 
Washburn  (Reference  1)  the  probability  that  area  A  contains  n  targets  is 

P(NA=n)  n=0,l,2 . 

This  is  the  Poisson  distribution  and  the  field  of  targets  is  termed  a  "Pois¬ 
son  field." 

The  probability  that  exactly  n  targets  are  detected  given  that  n 
targets  are  present  is 


where 


Pn  is  the  probability  that  exactly  n  targets  are  detected 

W  Is  the  sweepwldth 

7  Is  the  effective  search  speed 

A  Is  the  area  to  be  searched 

T  Is  the  time  duration 

The  "effective  speed"  Is  the  average  relative  speed  between 
searcher  and  target  assuming  the  angle  between  their  velocity  vectors  are 
uniformly  distributed  between  0  and  2*. 
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By  changing  the  variable  of  Integration  to  if)  =  (  it  -  9 )/2  and 
Introducing  $1noa  2>fuv/ { u+v) ,  the  equation  for  V  reduces  to: 

V  =  |(u+v)E(sin  a) 

where  E(slno)  is  the  complete  elliptic  integral  of  the  second  kind  and  sin  a 
Is  as  defined  above. 

The  probability  that  there  are  no  undetected  targets  in  Area  A 
after  search  of  duration  T  becomes  the  summation  of  the  Individual  terms 
PnP(NA=fO  and  thus: 

P(no  undetected  targets  In  A  after  search  of  duration  T). 


letting  p  *  1-e"^-,  multiplying  each  term  by  e'PPAePP*  and 
regrouping  we  get: 


p  «  ^  (P-P*)n  e"^  e_yA{1'p) 

11=0 

OD 

,  e-PA{l-p)  ^  e-VV* 

n=o 

where  the  Infinite  sum  Is  over  the  individual  terms  of  the  Poisson 
distribution  and  Is  Identically  equal  to  1.  The  expression  therefore  reduces 
to: 
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p  =  e-pA(l-p) 
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This  algorithm  had  not  been  programmed  at  the  time  work  was 
halted.  In  order  to  bring  this  development  to  completion,  the  sweep  width 
must  be  computed  and  V  must  be  calculated.  Models  for  computation  of  these 
two  factors  are  available  at  Summit  Research  and  the  remainder  of  work  on 
this  algorithm  would  be  to  Integrate  all  the  components  to  allow  for  a  more 
highly  automated  system  eliminating  the  necessity  for  table  look-ups  and  hand 
calcul atlons. 

SSN  EFFECTIVENESS  IN  MULTIPLE  ENGAGEMENT  SCENARIO 

No  additional  work  past  concept  formulation  occurred  because  of  the 
termination  of  this  task.  Presented  below  Is  the  concept  formulation  as  pro¬ 
posed. 

For  multi -engagement  models  there  are  two  obvious  approaches  that 
could  be  followed.  One  approach  would  be  to  simulate  successive  engagements 
in  a  Monte  Carlo  simulation.  An  alternative  approach  Is  to  use  Markov  anal¬ 
ysis.  Monte  Carlo  simulations  are  time  consuming  and  expensive  both  to 
create  and  to  run.  Markov  processes,  where  adaptable,  are  simple  to  formu¬ 
late  and  to  obtain  results. 


The  most  significant  requirement  for  formulating  the  multi - 
engagement  situation  as  a  Markov  process  Is  the  Identification  of  independent 
states.  In  a  one-on-one  engagement.  Independent  states  would  reflect  reduced 
systems  capabilities  on  either  side,  or  a  change  In  Information  state.  The 
following  are  examples  of  states  for  a  submarine  versus  submarine  engagement: 

•  all  systems  functioning  and  unalerted,  for  both  sub¬ 
marines; 
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•  all  systems  functioning  and  alerted,  for  both  subma¬ 
rines; 

•  one  submarine  has  reduced  capability  due  to  damage, 
other  fully  operational; 

•  both  submarines  have  reduced  capabilities  due  to  sys¬ 
tems  damage; 

•  etc. 

Terminal  states,  where  the  engagement  ends,  would  have  to  be  specified. 
These  states  would  reflect  the  situations  In  which  both  parties  no  longer 
want  to  pursue  the  engagement  any  further.  These  would  depend  on  the  objec¬ 
tives  of  both  sides.  For  some  scenarios,  damage  to  own  ship  or  to  the  oppo¬ 
nent  will  be  sufficient  for  disengagement.  In  other  situations,  an  engage¬ 
ment  may  continue  until  one  submarine  Is  disabled  or  destroyed,  or  the 
objective  is  accomplished. 

Once  the  states  are  specified  it  is  necessary  to  determine  the 
probability  of  ending  In  that  state  as  an  outcome  of  an  engagement.  Each  set 
of  transition  probabilities  must  be  calculated  assuming  the  engagement  starts 
with  the  combatants  In  each  of  the  initial  states.  Therefore,  the  engagement 
analysis  must  be  performed  for  each  of  the  states. 

While  a  good  deal  of  analysis  must  be  performed  for  this  approach, 
in  general.  It  is  easier  and  faster  than  using  Monte  Carlo  simulations.  The 
key  to  a  successful  use  of  Markov  techniques  Is  the  Identification  of 
states.  Each  state  must  reflect  a  significant  change  of  capabilities.  If 
too  many  states  are  specified,  the  analysis  and  computation  can  become  over¬ 
whelming.  If  too  few  are  selected,  the  sensitivity  of  results  can  be  lost. 
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This  same  approach  can  be  used  for  multiple  unit  scenarios.  This 
approach  was  successfully  used  In  the  FORCE  MIX  study.  The  states  for  a 
multiple  unit  engagement  reflect  decreases  In  the  number  of  units  available. 

The  more  complex  the  problem,  the  greater  the  advantage  of  Markov 
analysis  over  Monte  Carlo  simulations. 
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APPENDIX  A 


Presented  In  this  appendix  is  the  indexing  scheme  for  reconstruct¬ 
ing  the  2Np(N+l)  possible  tracks  used  In  the  barrier  model.  Figure  A-l 
contains  the  FORTRAN  code  pertinent  to  this  reconstruction.  This  code 
assumes  the  existence  of  a  three  dimensional  array,  DQS(L,KM,II).  An  element 
of  this  array  Is  the  probability  of  not  detecting  for  the  segment  of  relative 
track  ending  at  grid  point  (L,KM).  The  II  Index  determines  whether  this  seg¬ 
ment  enters  the  grid  point  from  the  right  or  left.  That  is,  if  11*1  then  the 
segment  Is  from  (L+l.KM+l)  to  (L,KM)  and  If  11=2  It  Is  from  (L+1,KM-1)  to 
(L.KM). 

From  Figure  A-l,  lines  1-3  set  the  indices  for  the  loops  (lines 
4,5)  for  all  relative  starting  positions  of  barrier  and  transitor  subma¬ 
rines.  There  are  (NSEG+1)  starting  positions  for  the  transitor  submarine 
(line  4)  across  the  barrier  width  and  (NPSEG+1)  starting  position  for  the 
barrier  submarine  (line  5).  Indices  In  lines  2  and  3  center  the  patrol  seg¬ 
ments  over  the  barrier  width. 


Code  In  lines  6-8  determines  the  Index  value  for  line  10.  This 
Index  (JU)  Is  associated  with  the  possible  directions  of  travel  of  the  bar¬ 
rier  submarine  at  Its  Initial  starting  position.  The  Index  OU  will  equal  2 
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1.  IB  =  NSEG+1 

2.  HI  *  (NSEG-NPSEG)/2.0  +1 

3.  M2  -  Ml  +  NPSEG 

4  DO  500  K  -  1,  IB 

5.  DO  500  M  =  Ml,  M2 

6.  JU  =  2 

7.  IF(M.EQ.Ml)  JU  =  1 

8.  IF(M.EQ.M2)  JU  =  1 

9.  1*1 

10.  00  600  J  =  1,  JU 

11.  IC  =  IC  +  1 

12.  TQ  *  1.0 

13.  MJ  *  M 

14.  DO  700  L  -  1,  LB 

15.  II  *  1 

16.  IF(MJ .EQ.M2)  I  *  -1 

17.  IF(MJ.EQ.Ml)  I  *  1 

18.  KM  =  MJ  -K  +  IB 

19.  IF(I.EQ.-l)  II  =  2 

20.  DQSS  *  DQSa.KM.II) 

21.  TQ  -  TQ*DQSS 

22.  MJ  *  MJ  +  I 

23.  700  CONT. 

24.  TBEFF  *  1.0  -  TQ 

25.  BEFF  -  BEFF  +  TBEFF 

26.  I  -  -1 

27.  600  CONT 

28.  500  CONT 

29.  BEFF  -  BEFF/IC 


Figure  A-l  CODE  FOR  RECONSTRUCTING  RELATIVE  TRACKS 
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unless  the  barrier  submarine  starts  at  one  of  the  endpoints  of  its  patrol 
width  (Ml  or  M2)  and  then  only  one  direction  Is  possible. 

The  Index  I  determines  the  relative  direction  of  thtf  barrier  sub¬ 
marine  for  each  segment  of  track.  It  is  +1  when  the  barrier  submarine  is 
moving  to  the  right  and  -1  when  moving  to  the  left.  The  value  of  I  deter¬ 
mines  the  index  II  in  the  DQS  array  mentioned  previously.  The  first  time 
through  the  J  loop  (lines  10-27)  the  value  is  set  to  1,  (line  9)  the  second 
time  through  the  value  Is  set  to  -1  (line  26)  for  the  different  Initial 
directions  of  travel  of  the  barrier  submarine.  If  the  initial  barrier  sub¬ 
marine  position  is  at  the  extreme  right  of  its  patrol  then  only  an  Initial 
position  applies  (I  =  -1)  and  line  16  correctly  sets  I  -  1  (line  9)  to  a  -1 
for  this  track. 

The  total  number  of  tracks  are  accumulated  by  the  expression  in 
line  11  which  was  initialized  to  zero  before  starting  this  reconstruction. 
Line  12  initializes  the  value  of  the  probability  of  not  detecting  along  any 
relative  track.  Line  13  establishes  a  dummy  variable  that  can  modified  with¬ 
out  changing  the  value  of  the  M  index  from  line  5. 

The  probability  of  not  detecting  along  one  of  the  relative  tracks 
Is  computed  by  lines  14-23.  As  the  transitor  moves  from  its  initial  to  its 
final  position  it  traveres  LB  range  Increments.  Depending  on  the  initial 
position  of  the  transitor  across  the  barrier  width  and  the  initial  position 
and  direction  of  the  barrier  submarine,  (initial  values  of  K,  M,  and  I).  The 
track  is  reconstructed  by  varying  the  KM  index  (line  18)  by  *1  for  each  range 
Increment  In  L  depending  on  whether  the  barrier  submarine  would  be  moving  to 
the  left  or  right  during  this  Increment.  Lines  16  and  17  determine  whether 
the  endpoints  of  the  patrol  are  reached  and  a  change  of  direction  Is  neces¬ 
sary.  Line  21  Is  the  probability  of  not  detecting  along  the  track  (P(not  de¬ 
tecting)  *  jilbql)  • 
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Line  24  Is  the  probability  of  detecting  along  the  track  and  line  25 
accumulates  this  probability,  which  Is  eventually  divided  by  the  total  number 
of  tracks  IC  (line  29)  yielding  the  mean  probability  of  detection  for  the 
moving  barrier. 
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